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ABSTRACT 

We explore the question of whether low and moderate luminosity Active Galactic Nuclei (AGNs) 
are preferentially found in galaxies that are undergoing a transition from active star formation to 
quiescence. This notion has been suggested by studies of the UV-optical colors of AGN hosts, which 
find them to be common among galaxies in the so-called 'Green Valley', a region of galaxy color space 
believed to be composed mostly of galaxies undergoing star-formation quenching. Combining the 
deepest current X-ray and Herschel. PACS far-infrared (FIR) observations of the two Chandra Deep 
Fields (CDFs) with redshifts, stellar masses and rest-frame photometry derived from the extensive and 
uniform multi-wavelength data in these fields, we compare the rest-frame U — V color distributions 
and SFR distributions of AGNs and carefully constructed samples of inactive control galaxies. The 
UV-to-optical colors of AGNs are consistent with equally massive inactive galaxies at redshifts out to 
z ~ 2, but we show that such colors are poor tracers of star formation. While the FIR distributions 
of both star-forming AGNs and star-forming inactive galaxies are statistically similar, we show that 
AGNs are preferentially found in star-forming host galaxies, or, in other words, AGNs are less likely 
to be found in weakly star-forming or quenched galaxies. We postulate that, among X-ray selected 
AGNs of low and moderate accretion luminosities, the supply of cold gas primarily determines the 
accretion rate distribution of the nuclear black holes. 



1. INTRODUCTION 

Accreting black holes are arguably the most efficient 
engines of energy production in the Universe. The deep 
gravitational wells of supermassive black holes (SMBHs) 
allow the extraction of ~ 10% of the rest mass-energy of 
the material that falls into their horizons, which, through 
accretion processes, is ultimately converted into the elec- 
tromagnetic and mechanical output that power Active 
Galactic Nuclei (AGNs). This energy can, in turn, es- 
cape into the environs of the AGN host galaxy, affecting 
material on large scales. Such AGN 'feedback' has many 
potential effects on galaxy physics and evolution: reg- 
ulation of the circum-nuclear environment and galactic 
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star- formation (SF), gaseous outfiows, the distribution 
of metals, enrichment and heating of gas in the circum- 
galactic medium. There is much observational evidence 
for the direct effects of SMBH activity on their host 
galaxies, either in the form of AGN-driven outfiows (Holt 
et al. 2008; Alexander et al. 2010; Feruglio et al. 2010; 
Fischer et al. 2010; Greene et al. 2011; Rupke & Veilleux 
2011; Sturm et al. 2011; Cano-Dfaz et al. 2012; Harrison 
et al. 2012b) or bubbles blown in the host atmospheres of 
massive galaxy haloes by powerful radio jets (McNamara 
& Nulsen 2007; Gitti et al. 2012). 

Proposed evidence for the widespread action of SMBH 
feedback on the star- formation histories of galaxies comes 
from suggestions that most AGNs occupy a preferred 
population of host galaxies, those that are undergoing 
a transformation from active steady star-formation to a 
final state of quiescence (e.g., Nandra et al. 2007; Martin 
et al. 2007; Schawinski et al. 2010). These two phases 
form the basis of the well known color-bimodality of 
galaxies - most galaxies from the local Universe to z ^ 3 
lie in two relatively distinct parts of the color-magnitude 
or color-mass diagram, the Blue Cloud and the Red Se- 
quence. Galaxies that lie at intermediate colors on these 
diagrams are believed to be transitioning between the 
two populations, through a region known as the 'Green 
Valley' (e.g., Faber et al. 2007). The preponderance of 
Green Valley galaxies among AGN hosts has been taken 
as evidence that low-to-moderate luminosity AGNs are 
responsible for the quenching of star-formation in galax- 
ies, mediated through feedback from the SMBH. The 
closely associated morphological evolution of transform- 
ing galaxies (Driver et al. 2006; Kauffmann et al. 2006; 
Franx et al. 2008; Cheung et al. 2012) is taken to be evi- 
dence that a substantial fraction of AGNs are related to, 
and perhaps triggered by, galaxy mergers which may be 
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responsible for the formation of 'Red and Dead' ellipti- 
cals. 

In recent work, the notion that AGN prefer quench- 
ing hosts has come under greater scrutiny. The impor- 
tance of stellar mass selection effects for the interpre- 
tation of AGN colors seems to suggest that AGN host 
colors are not radically different from similarly massive 
inactive galaxies (Xue et al. 2010; Rosario et al. 2013). A 
previous study by our team (Santini et al. 2012b) found 
that mean SFRs are enhanced in X-ray AGN over in- 
active galaxies of the same stellar mass, with tentative 
evidence that the enhancement was caused by a lower 
fraction of quiescent galaxies among AGN hosts, rather 
than a boost in the SFRs of star-forming AGNs. How- 
ever, given the relatively shallow far-infrared data used 
in this earlier work, this notion could not be tested ex- 
tensively. 

In this paper, we critically examine the evidence that 
AGN arc preferentially in galaxies that are quenching or 
undergoing a slow-down of their global star-formation. 
Rather than relying only on optical or UV tracers of 
star- formation such as the color-mass diagram (CMD), 
we employ instead the most sensitive far- Infrared (FIR) 
data currently available from the Herschel Space Tele- 
scope, together with a large sample of AGNs selected 
from the deepest cxtragalactic X-ray surveys on the sky 
- the two Chandra Deep Fields. In §2, we present the sur- 
veys, selection and datasets. In §3, we compare the use 
of FIR tracers and UV-optical colors in studies of SF. In 
§4 and 5, we re-examine the evidence that AGNs lie pref- 
erentially in the Green Valley, then fold in information 
about the FIR luminosities and detection rates of AGNs 
toward exploring the question of whether AGNs are in 
quenching galaxies. Our results are discussed in §6. In 
this work, we assume a A-CDM concordance cosmology 
with Qa = 0.7, flrn = 0.3 and Hq = 70 km s'^Mpc-^ 

2. DATASETS AND SAMPLE SELECTION 

2.1. Sample Selection and Datasets 

Gospatial with the two GOODS survey fields (Gi- 
avalisco et al. 2004), the Chandra Deep Fields (CDFs) 
are the deepest pencil-beam X-ray surveys in the sky. In 
GOODS-North, the CDF-North (CDF-N) X-ray catalog 
comprises 503 sources from a total exposure of 2 Msec 
(Alexander et al. 2003), while in GOODS-South, the new 
4 Msec CDF-South (CDF-S) X-ray catalog consists of 
740 sources (Xue et al. 2011). We have extensively char- 
acterized the data and catalogs in both fields, in which 
careful associations have been made with optical and 
near-IR counterparts, using, where possible, probabilistic 
crossmatching models (Luo et al. 2010; Xue et al. 2011). 
In addition to the dec^p X-ray data, the wealth of deep 
spectroscopy and irnilti-wavelength photometric data in 
the GOODS fields have enabled accurate spectroscopic or 
AGN-optimised photometric redshifts to be determined 
for the majority of the X-ray sources (e.g., Szokoly et al. 
2004; Luo et al. 2010). We estimate absorption-corrected 
hard-band X-ray luminosities {Lx) of sources with red- 
shifts using spectral modeling techniques (Bauer et al. 
2004) . As a result of the small area and great depth of 
the CDF exposures, most X-ray sources are low or mod- 
erate luminosity AGNs - only ~ 5% of the sources have 
logXx(2-10 keV) > 44 erg s~^. These have equivalent 



AGN bolometric luminosities to the Seyfert galaxy popu- 
lation found in the local Universe. In this work, we only 
consider sources with logL^ > 42 erg s~^, to prevent 
contamination from powerful starbursts, in which emis- 
sion from X-ray binaries can potentially overpower the 
emission from nuclear activity in such faint systems. 

We employ multiwavelength galaxy catalogs for the 
two GOODS fields to define a general galaxy sam- 
ple, the properties of which we will compare to the 
AGNs. In GOODS-S, we use the updated GOODS- 
MUSIC database (Santini et al. 2009; Grazian et al. 
2006), while in GOODS-N we use a catalog developed 
for the PEP team using similar methodologies (Berta 
et al. 2010, 2011) The former catalog selects galax- 
ies with observed magnitudes in the HST F850LP band 

< 26 or in the ISAAC band < 23.5, while the lat- 
ter is primarily selected to have K < 24.2. In order to 
exclude a surfeit of faint sources with inaccurately red 
colors and masses, we apply an additional cut of F850LP 

< 26 in the GOODS-N catalog. For galaxies with no cur- 
rent spectroscopic redshifts, photometric redshifts were 
determined by fitting multiwavelength photometry us- 
ing PEGASE 2.0 templates (Fioc & Rocca-Volmerange 
1997) in GOODS-S or using the EAZY code (Brammer 
et al. 2008) in GOODS-N. For details on the catalog 
preparation, characterization and photometric redshift 
estimation, we refer the reader to Santini et al. (2009) 
and Berta et al. (2010) for GOODS-S and GOODS-N 
respectively. 

While AGN are selected by their X-ray emission, we de- 
fine our 'inactive' galaxy population as all galaxies that 
are undetected in X-rays (excluding even those which 
have \ogLx < 42 erg s~^) and with mid-IR (MIR) col- 
ors (i.e, based on Spitzer/IRAC photometry) that do 
not satisfy the AGN selection criteria of Donley et al. 
(2012). In practice, only a very small fraction of the 
general galaxy population are rejected on the basis of 
these criteria. These rejected objects, however, tend to 
be in massive galaxies and could potentially sway the 
statistics of SF comparisons among such systems by an 
inordinate degree. 

We have developed a custom technique for the esti- 
mation of stellar masses (M*) in AGNs by linearly com- 
bining galaxy population synthesis model templates and 
AGN SED templates to fit multiwavelength photometry. 
For inactive galaxies, we perform a minimization of 
Bruzual & Chariot (2003) synthetic models, assuming a 
Salpeter Initial Mass Function (IMF) and parameteriz- 
ing the star formation histories as exponentially declining 
laws. For AGNs, we also include an AGN template from 
Silva et al. (2004), which accovmts for a variable fraction 
of the total light of the galaxy. The AGN template re- 
flects the classification of the X-ray source, derived from 
information about its SED and spectrum, where avail- 
able. For sources classified as Type I (broad lines in the 
spectrum, clear AGN contribution in the rest-frame opti- 
cal and UV), a Seyfert 1 SED was used, while for the rest, 
a Seyfert 2 template was used if the estimated X-ray ab- 
sorption column Nh < 10^* cm^^, and a Compton-thick 
template for more heavily absorbed systems. For further 
details, performance evaluations and tests of the method. 

Available on the PEP public release page. 
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Fig. 1. — Color-mass Diagrams (CMDs) of rest-frame U — V color against stellar mass (M*). All X-ray sources and IRAC-selected AGNs 
have been excluded from these diagrams. The general galaxy population taken from the full galaxy catalogs in both GOODS fields are 
shown as small black points, with the distribution oi U — V color for galaxies with M«> 1O^'^M0 (dotted line) shown as a black histogram 
on the right axis of each panel. A color bimodality is seen at all redshifts. The location of the 'Green Valley' is shown by solid green lines 
which mark the locus of minimum galaxy density in this diagram (see §3 and Table 1). FIR-detected galaxies (for a definition, see §2) are 
shown as large colored points and their U — V distribution is shown on the right axis of each panel as a colored histogram. Both histograms 
are normalized to have the same peak value. Despite lying on or above the SF Mass Sequence, most FIR detected galaxies exhibit U — V 
colors that are intermediate between the Red and Blue sequences, i.e, they lie in the Green Valley and show a weak color bimodality. This 
is primarily because of their high stellar masses. 

S and 0.93/2.04 mJy at 100/160 fim in GOODS-N. The 
GOODS-S maps are « 80% deeper than the GOODS-N 
maps and probe further down the FIR luminosity func- 
tion at all redshifts (Magnelli et al. 2013). 

For practical purposes, we use the monochromatic lu- 
minosity of a galaxy at 60 /J,m rest (Leo) as a measure 
of its FIR luminosity. The PACS bands cover this rest- 
frame wavelength over much of the redshift range probed 
in this work and we estimate Lgo from a simple log- 
linear interpolation of PACS measurements in bands that 
bracket 60 /im in the rest-frame. The use of Lgo obvi- 
ates the need to apply an uncertain correction between 
monochromatic and total FIR luminosities. Neverthe- 
less, in order to relate Lgo to the properties of the pop- 
ulation of star-forming galaxies from the existing litera- 
ture, such as the SF 'Mass Sequence' or 'Main Sequence' 
(MS), we adopt the following relationship between Lgo 
and SFR: 



TABLE 1 
Parameters of the Green 
Valley (GV) lines in the 
Color-Mass Diagram 



Redshift interval 


UVg 


a 


0.5-1.0 


0.60 


0.17 


1.0-1.5 


0.50 


0.15 


1.5-2.0 


0.45 


0.15 



Note. — UVg and a are de- 
fined such that (U — V)gv = 
alog(M*,g) -I- C/Vg, where {U - 
V)gv is the color of the GV line 
and Mt^g is the stellar mass in 
units of 10^ solar masses. 



we refer the reader to Santini et al. (2012b). 

Our far-infrared data are composed of maps at 70 /im, 
100 /.tm and 160 /.tm from a combination of two large Her- 
schel/PACS programs: the PACS Evolutionary Probe 
(PEP), a guaranteed time program (Lutz et al. 2011) and 
the GOODS-Herschel key program (Elbaz et al. 2011). 
The combined PEP-hGH (PEP/GOODS-Herschel) re- 
ductions are described in detail in Magnelli et al. (2013). 
While data at 100 and 160 /im are available in both 
fields, an additional deep map at 70 /im is also avail- 
able in GOODS-S. The PACS 160, 100 and 70 //m fluxes 
were extracted using sources from archival deep Spitzer 
MIPS 24 /xm catalogs as priors, following the method 
described in Magnelli et al. (2009); see also Lutz et al. 
(2011) for more details. 3cr depths are 0.90/0.54/1.29 
mJy at 70/100/160 /zm in the central region of GOODS- 



LiR — fcEoi X Lgo 
SFR^l.72 X 10-^° 



in. 



(1) 
(2) 



fc Eoi is the conversion factor between Lgo and the to- 
tal IR luminosity integrated over 8-1000/fm [Liji) from 
Chary & Elbaz (2001). This factor is IR luminosity de- 
pendent, but only varies slightly over the range 1.6-2.5. 
The conversion factor between SFR and L/^ is taken 
from Kennicutt (1998), which also assumes a Saltpeter 
IMF. 
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Fig. 2. — Specfic Star Formation Rate (SSFR) against stellar mass (M,) of PACS-detected inactive galaxies. Points are colored by the 
A{U — V), the U — V offset of the galaxy from the Green Valley (GV), defined by straight lines in the CMD (§3). Galaxies that lie in the 
GV have A{U — V) close to zero and are colored as green points in the Figure. PACS-detected galaxies in the GV lie on or around the 
Mass Sequence (shown as a dashed lines) and are typically massive, while those in the Blue Cloud are typically at lower masses and lie well 
above the Mass Sequence. Galaxies in the GV correspond to normal massive SF galaxies at all redshifts in this study. 



3. THE COLORS AND MASSES OF FIR-DETECTED 
INACTIVE GALAXIES 

We begin by first examining the relationship between 
the optical color of galaxies and the star-formation rate 
or star-formation luminosity, as traced by FIR emission. 
For this, we turn to the color-mass diagram (CMD) as a 
diagnostic tool. In Fig. 1, we plot the rest- frame U — V 
color against the stellar mass of galaxies, separately 
for the two GOODS fields. Consistent with several pre- 
vious studies in these and other extragalactic fields (e.g., 
WiUmer et al. 2006; Wyder et al. 2007; Taylor et al. 2009; 
Brammer et al. 2009; Whitaker et al. 2011), galaxies tend 
to cluster in two well-defined regions of the diagram: the 
Red Sequence and the Blue Cloud. These structures form 
the basis of the well known color-bimodality of galaxies, 
which has been revealed to z > 2 (Brammer et al. 2009). 
In such optical- or NIR-selected photometric catalogs, 
the mass limit, revealed by the sharp boundary in the 
density of galaxies at the lower mass end of the diagram, 
is redshift and color dependent. At z > 1.5, the mass 
limit is high enough among red galaxies that it blurs the 
definition of the color bimodality. In addition, at these 
redshifts, a population of extremely dust-reddened galax- 
ies is also seen, leading to a tail of very red colors among 
high mass galaxies. 

In between the dense Red Sequence and Blue Cloud, 
galaxies have intermediate colors, which has led to the 
popular name for this area of the CMD: the Green Valley 
(GV). Taking the GV as a minimum in the density dis- 
tribution of galaxies on the CMD at a given stellar mass, 
we construct sloped lines on the CMD which separate 
the Red Sequence from the Blue Cloud and define the 
location of the GV for all our subsequent analysis. The 
slope and normalization of the lines were determined by 
eye to yield the most well-defined separation between the 
Red Sequence and Blue Cloud in each redshift bin. We 
tabulate the GV lines for each bin in Table 1 and plot 
them in Fig. 1 using solid green lines. 

It has been suggested, based on the low density of the 
GV, that most galaxies here are going through a rela- 



tively rapid phase of SF quenching, resulting in a net 
fiow of galaxies from the SF Blue Cloud to the quies- 
cent Red Sequence (Faber et al. 2007; Martin et al. 2007; 
Brammer et al. 2011). In other words, GV galaxies are 
believed to have lower specific SFRs compared to nor- 
mal SF galaxies, which are defined to lie on the SF Mass 
Sequence (Noeske et al. 2007; Elbaz et al. 2007; Santini 
et al. 2009; Rodighiero et al. 2011; Wuyts et al. 2011; 
Whitaker et al. 2012). 

Fig. 1 suggests otherwise. We have plotted in this Fig- 
ure, using large colored symbols, the locations of PACS 
detected inactive galaxies from the PEP-t-GH catalogs. 
In general, these galaxies lie on or above the MS at 
z > 0.5 and tend to be rather massive (Mh.> 10^" M©). 
It is clear from the plot that a substantial, number of 
FIR-bright galaxies lie in the GV. In fact, the main de- 
terminant of whether or not a galaxy is detected in PACS 
is its stellar mass, not its color. 

This can be examined in a different way using Fig. 2, 
where we have plotted the FIR-derived specific SFR 
(SSFR) of PACS-detected inactive galaxies against M*. 
The ridgeline of the SF MS at the central redshift of each 
redshift bin is shown in these plots as dashed lines, as de- 
termined recently by Whitaker et al. (2012). The points 
here are colored by A{U — V), the U — V color offset 
of the galaxy from the GV in the CMD, i.e the vertical 
offset of the galaxy from the green GV lines in Fig. 1. 

The steeper slope of the PACS-detected galaxies com- 
pared to the MS lines is due to the flux limit of the PACS 
photometry, which translates into a mass-dependent 
limit in SSFR at a given redshift. At lower M,, only 
galaxies that lie at progressively higher SFR above the 
MS are detectable in PEP-I-GH. This strongly affects the 
interpretation of the slope of the MS purely from FIR 
data, which is why we adopt MS relationships deter- 
mined, in this case, from a uniform study of deeper UV 
and 24/im photometry (Whitaker et al. 2012). We have 
changed the mass scaling in the relations of Whitaker 
et al. (2012), estimated using a Chabrier IMF, upwards 
by a factor of 1.74 to reflect our use of the Salpeter IMF 
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(Santini ct al. 2012a). 

At M^,> lO^"'^ Mq, wc sec that galaxies that lie around 
the MS have colors that are typical of the GV, i.e with 
A{U - V)!v (see also Whitaker et al. 2012; Salmi 
et al. 2012). In essence, this means that massive galax- 
ies with intermediate colors are not a special population 
that is quenching and moving away from a state of on- 
going SF, but, indeed, are the star-forming population 
at M*> IQiO-S M0. 

In later sections, we examine the PACS detection frac- 
tions of AGNs and inactive galaxies. It is worthwhile, 
at this stage, to note that PACS detected galaxies pri- 
marily lie on or above the MS, even though the MS is 
known to have a roughly symmetrical scatter of about 
0.3 dex in SSFR (Noeske et al. 2007; Rodighiero et al. 
2011). Concentrating only on PACS detected galaxies 
restricts us to galaxies that scatter above the MS, but 
one must remember that there are normal star-forming 
galaxies that nominally lie within the MS but are below 
our detection limit since they lie within the scatter below 
the MS. When discussing PACS detection statistics, we 
lump together quiescent and quenching galaxies, which 
are expected to lie well below the normal scatter of the 
MS, and 'weakly SF galaxies', those that lie in the low 
scatter of the MS. The reader should keep in mind that 
the population of PACS-undetected galaxies are not all 
quiescent or quenching, but will in addition include a 
sizable number of actively SF galaxies that are simply 
below our nominal PACS detection limits. 

4. THE COLORS OF X-RAY AGN: COMPARISON TO 
INACTIVE GALAXIES 

We have shown the rest-frame optical color is not an 
ideal tracer to identify a quenching population because 
its use as a measure of the SF properties of a galaxy de- 
pends on the stellar mass of the galaxy. Nevertheless, 
several previous studies have taken the location of AGNs 
in the CMD as evidence for a link between AGN activity 
and the transformation of galaxies (e.g., Nandra et al. 
2007; Martin et al. 2007; Schawinski et al. 2010). The 
reason for this can be seen in Fig. 3. The AGNs, shown 
as large colored points, scatter mostly about the GV lines 
that we defined in §3, with a small scatter to lower masses 
and colors in the Blue Cloud. In many ways, the approx- 
imate distribution of AGN in the CMD mirrors that of 
PACS-detected inactive galaxies. However, this is almost 
completely because both AGNs and FIR-bright galaxies 
are typically quite massive. 

Indeed, it is this strong tendency for X-ray detected 
AGNs to lie in massive host galaxies that mostly deter- 
mines the location of the AGNs in the CMD. Recent 
studies suggest that, once the particular mass distribu- 
tion of AGN hosts is taken into account, the colors of 
AGNs and inactive galaxies are very similar (Silverman 
et al. 2009; Xue et al. 2010; Cardamone et al. 2010; 
Rosario et al. 2013). The small sizes of AGN samples, 
generally only a few % of galaxies, make it difficult to 
identify statistically robust differences in the colors of 
AGNs and inactive galaxies. Here we develop a method 
that builds on the much larger sample of inactive galax- 
ies to test the following null hypothesis: AGN hosts are 
drawn randomly from the population of massive galaxies 
and share the SFRs and colors of the parent sample. 

For each AGN, we choose, at random and allowing 



duplicates, an inactive galaxy in the same redshift bin 
and with a stellar mass within ±0.1 dex of the mass of 
the AGN host galaxy. In this way, we arrive at a sam- 
ple of inactive galaxies which are equal in number and 
mass distribution as the AGNs in the redshift bin. The 
mass tolerance is smaller than the typical error in stellar 
masses from our SED fits, ensuring essentially identical 
mass distributions of AGNs and inactive galaxies. From 
this comparison sample, we derive a distribution in color. 
We repeat this process for a total of 1000 trials. This 
bootstrap approach gives us the typical uncertainty in 
the color distribution for inactive galaxies that share the 
mass distribution of AGNs. 

In Fig. 4, we compare A{U — V), the U — V color off- 
set from the Green Valley, of AGNs and inactive galaxies. 
The distribution of A{U — V) for the AGNs are shown 
as a open histogram, plotted over the distributions of an 
equal number of inactive galaxies with the same redshift 
and stellar mass range as the AGNs. Using our bootstrap 
procedure, we obtain la and 2a uncertainties on the in- 
active galaxy distributions, shown as dark and light grey 
shaded regions in the Figure. The colors of AGNs are 
statistically different from those of the inactive galaxies 
only if the open histogram strongly or consistently devi- 
ates from the darkly shaded regions over a few or more 
bins in A{U — V). Note that systematic differences in the 
distributions due to cosmic variance are not represented 
in the uncertainties of these shaded histograms, though 
they can be important in small fields like GOODS and 
among the massive galaxies considered here. A compar- 
ison using both fields helps in this regard. 

A quick examination shows that, in almost all redshift 
bins, there are very minor differences between AGNs and 
inactive galaxies. At the location of the Green Valley, the 
AGNs are slightly more common than inactive galaxies, 
at the level of l-1.5a in five of six panels, but this differ- 
ence is quite small and may be attributable to enhanced 
line emission or recombination continuum emission in 
AGNs, since nuclear activity can frequently produce ex- 
tended, highly ionized, emission line regions. At about 
the same low level of significance, there is an enhance- 
ment of inactive galaxies over AGNs in the Red Sequence 
{A{U — V) in the range 0.2 — 0.6 mag) . Perhaps the most 
significant difference between the two distributions is not 
near their peak, but is in the extreme blue wing, where 
AGNs are more common than inactive galaxies. This is 
likely due to the small fraction (a few %) of AGNs which 
contain bright blue nuclear point sources which can dom- 
inate the integrated light of the system. These will have 
very blue colors for their stellar mass, a result of nuclear 
contamination in the optical and UV bands. 

5. THE FIR-DERIVED SFRS OF X-RAY AGN: 
COMPARISON TO INACTIVE GALAXIES 

We now turn to a far better tracer of SF in galaxies: 
the FIR luminosity. Among massive galaxies, such as X- 
ray selected AGN hosts, the FIR luminosity, an excellent 
tracer of the total dust-reprocessed UV light from star- 
forming regions, accounts for almost all the SF in the 
galaxy the UV escape fractions in such galaxies is quite 
low and corrections for un-reprocessed SF luminosity are 
at the level of a few % (Pannella et al. 2009; Reddy et al. 
2010; Whitaker et al. 2012). 

The relative SF properties of AGN hosts and inactive 
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Fig. 3. — Color-mass Diagrams (CMDs) of rest-frame U — V color against stellar mass (M*) comparing AGNs and inactive galaxies. 
The general galaxy population of inactive galaxies are shown as small black points, with the distribution of U — V color for galaxies with 
M»> 1O®'^M0 (dotted line) shown as a black histogram on the right axis of each panel. X-ray selected AGNs are shown as large colored 
points and their U — V distribution is shown on the right axis of each panel as a colored histogram. Both histograms are normalized to 
have the same peak value. The location of the 'Green Valley' is shown by solid green lines as in Fig. 1. AGNs exhibit typically high stellar 
masses, similar to FIR-detected star-forming galaxies (compare this plot to Fig. 1). 



galaxies can be understood more accurately by compar- 
ing their FIR luminosity distributions. Studies of the 
FIR SEDs of AGNs show that the contribution to Lgo 
from dust emission heated by the active nucleus is low 
and generally negligible, except in a few cases of lumi- 
nous AGNs in weakly SF galaxies (Netzer et al. 2007; 
MuUaney et al. 2011; Rosario et al. 2012). We can di- 
rectly verify this for FIR-bright AGNs in the CDF fields, 
as done in Fig. 5, where we compare measured Lgo of 
AGNs detected in FAGS with the predicted FIR lumi- 
nosity derived from Lx based on model relationships of 
pure AGN-heated dust. These relationships, described 
in Rosario et al. (2012), have SED shapes covering the 
scatter found among local X-ray bright AGNs (MuUaney 
et al. 2011), while the nuclear X-ray luminosity is linked 
to the IR luminosity following the tight correlation of 
Gandhi et al. (2009). From the Figure, it is clear that 
the AGN luminosities of most of the sources in our sam- 
ple are too low to substantially affect the FIR luminosity. 

Another simple test is to compare the IR colors of 
AGNs and inactive galaxies; if the FIR colors of AGNs 
are significantly warmer than inactive galaxies (i.e, with 
higher relative flux in bluer bands), then AGN heating 
may be responsible for a non- negligible fraction of Lqq. 
In Fig. 6, we plot the MIPS 24 ^lm to PACS 160 /im 
color against the PACS 100 fj,m to 160 ^m color (both 
expressed as magnitudes) of PACS detected AGNs and 
inactive galaxies. This plot may be directly compared to 
Fig. 2 of Rosario et al. (2012), which studied the IR colors 
of sources in the COSMOS survey in similar fashion. The 
FIR color distributions of AGNs and inactive galaxies 
(histograms along the Y-axis in all panels) are formally 
indistinguishable, though, expectedly, AGN emission can 
influence the MIR, as evidenced by the scatter of AGN 



points to bluer 24-to-160 fim colors. For essentially all 
of our PACS-detected AGNs, we can safely assume that 
Lqo is a measure of the SFR. 

In Fig. 7, we compare ALgo, the offset of a galaxy from 
the Whitaker et al. (2012) Mass Sequence, for AGNs and 
mass-matched inactive galaxies. The known rcdshift and 
stellar mass of each object determines the corresponding 
MS value; then the Chary & Elbaz (2001) IR SED li- 
braries and the relationship from Kennicutt (1998) are 
used to connect the MS to Lqq (§2). In practice, the 
choice of comparing simple Lqq distributions or ALgg 
distributions does not alter our basic result; we use the 
offset from the MS simply to scale out any stellar-mass 
dependent trend which would otherwise broaden the dis- 
tributions. As for Fig 4, we use a bootstrap sampling 
procedure, randomly selecting a set of inactive galaxies 
matched to the AGNs in M*, repeated a thousand times, 
to arrive at statistical uncertainties in the Lgg distribu- 
tions of the inactive galaxies. Dark/light grey regions 
show Icr and 2a uncertainties respectively. 

Unlike in A{U — V), there is a great deal of scatter in 
the inactive galaxy population in terms of FIR luminos- 
ity, as evinced by 2a uncertainties, which can span from 
zero to almost twice the number of AGNs at some values 
of ALgQ. Despite this, we find that the AGN distribu- 
tions are generally consistent with the inactive galaxy 
distributions at the 1-1. Scr level, i.e, the solid histograms 
lie within or close to the dark-grey regions. There is a 
possible minor tendency for an excess of PACS-detected 
SF galaxies at or below the ALqq— line compared to 
AGNs in the GOODS-S panels at z > 1. However, this 
is not seen in GOODS-N and may just be due to cos- 
mic variance. The broad consistency between the dis- 
tributions implies that any differences in SF properties 
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Fig. 4. — A comparison of the Green Valley eolor ofTset {A(U — V)) of X-ray selected AGNs and mass-matched inactive galaxies. The 
open histograms drawn with a solid black line show the distribution of A(U — V) for the AGNs, while shaded histograms correspond to the 
inactive galaxies. The statistical uncertainty in the distributions of the inactive galaxies are shown by the shading in the histograms - dark 
grey sections show the la uncertainty, due to the scatter in the population as well as small number statistics, while the light grey sections 
show 2a. The dashed line at A{U — V)= is the location of the center of the Green Valley. The AGNs show very similar distributions to 
ina<;tive galaxies, especially among the bluer stax-forming population, but are under-represented in the Red Sequence. 
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Fig. 5. — Leo plotted against the X-ray luminosity of the AGNs 
in both Chandra Deep Fields over the redshift range 0.5 < z < 
2.0. Different symbols are used to represent sources from the two 
different survey fields. The region between the two solid lines is 
where one may expect to find AGN with negligible star-formation 
and FIR SEDs dominated by dust heated by the nucleus (see §5 
and (MuUaney et al. 2011)). Essentially, all the AGN have FIR 
luminosities well above the pure-AGN region. It is highly unlikely 
that AGN emission powers the FIR luminosity of these sources, 
even among those at the luminous end. Note, the apparent wealc 
correlation between L2-10 and Leo is a consequence of Eddington 
bias and is not a real trend. 

between AGN hosts and inactive galaxies that lie on or 
above the SF Mass sequence are minor or non-existent, 
a result that is completely consistent with earlier work 
on the GDF X-ray selected AGN based on independent 



PEP and GOODS-Herschel data (Mullaney et al. 2012b; 
Santini et al. 2012b). Star-forming AGN hosts are drawn 
from the general population of SF galaxies. 

Galaxies with the greatest SF offset from the MS are 
known to be strong starbursts and frequently show highly 
disturbed morphologies consistent with being galaxy 
mergers (Wuyts et al. 2011). The Aigo distribution is 
dominated by galaxies close to the ridgcline of the MS 
and we have rather poor statistics for AGN and galax- 
ies at the highest Aigo 10), especially considering 
the substantial scatter shown in the distributions of in- 
active galaxies. Nevertheless, it is clear that AGNs are 
not strongly enhanced among starbursting galaxies, at 
least among the range of nuclear luminosities probed in 
the work. 

There are, however, some clear differences between the 
populations of AGN hosts and inactive galaxies which 
come to light once we consider sources that are not de- 
tected in the PACS maps. These sources typically lie 
below the MS and include quenched galaxies and weakly 
star-forming galaxies (i.e, normal SF galaxies that lie 
within the lower scatter of the MS). To understand this 
population, we compute the FIR 'non-detection' fraction, 
fnd which is the fraction of galaxies in each redshift bin 
that have no flux in both 100 and 160 /xm PACS maps 
at the 3a detection threshold of the catalogs, fnd is a 
measure of the relative number of weakly star-forming 
or quiescent galaxies in a population. 

In Fig. 8, we plot the histogram of fnd for inac- 
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Fig. 6. — Infrared colors (in magnitudes) of PACS-detected AGNs (red star points) and inactive galaxies (small black points) from both 
GOODS fields combined. 24 - 160 fim (a MIR-FIR color) is plotted on the X-axis and 100 - 160 /im (a FIR-FIR color) is plotted on the 
Y-axis. Distributions of 100 - 160 fim colors are shown as histograms on the Y-axes, with red/black colors for AGNs/inactive galaxies 
respectively. Both histograms are normalized to the same peak value. There is a small scatter of AGNs to low, blue MIR-FIR colors, due to 
the influence of AGN-heated dust on the 24 /im emission, but most AGNs lie in the area of the diagram occupied by inactive SF galaxies, 
implying that their FIR emission, if not the entire MIR-to-FIR SED, is dominated by SF- heated dust emission. 



tive galaxies determined from the 1000 realizations of 
the mass-matching bootstrap procedure described above. 
The median non-detection fractions for this population 
are approximately 80%. This fraction is roughly inde- 
pendent of redshift despite the evolving luminosity limit 
of the PACS data, which reflects the fact that galaxies 
at z ~ 2 are typically more star-forming, and hence FIR 
luminous, than at z ~ 0.5. In these insets, we show the 
value of fnd for the AGNs in that redshift bin, a single 
number, using a thick downward-pointing arrow. It is 
immediately clear that the AGNs are always, at all red- 
shifts and consistently in both fields, significantly more 
likely to be detected in PACS than the inactive galaxies. 
In other words, the fraction of quenching or quiescent 
galaxies hosting AGNs is considerably lower than simi- 
larly massive inactive galaxies, or, put differently, AGN 
are more likely to be in star-forming galaxies around or 
above the MS. This result is in stark contrast to the no- 
tion that AGN are preferentially in quenching systems. 
Controlling for the particular mass distribution of X-ray 
selected AGN hosts, we show that, in fact, AGN are 
more likely to be found in a galaxy that is forming stars, 
compared to one that is turning quiescent. This result 
strongly confirms the result from Santini et al. (2012b) 
that AGN have higher FIR detection rates, which leads 
to the enhancement in mean SFR among AGNs com- 
pared to normal galaxies. 

5.1. Biases and tests 

We have undertaken a series of tests to verify that this 
result is not due to possible biases in the data. Since 
the CDF sensitivity has a radial variation, the density of 
X-ray sources is higher in the central area of the field, 
while the comparison sample of background galaxies has 
a more uniform spatial distribution. Therefore, a pos- 
sible bias may arise if the equivalent PACS coverage, 
which varies across the PEP-I-GH maps, differs between 
the AGNs and inactive sample. An examination of the 
PACS data coverage of AGNs and inactive galaxies shows 



that they are quite consistent; this is because the parts 
of the GOODS fields with strong PACS coverage gradi- 
ents are restricted to the edges. We then repeated our 
analysis of FIR properties, while matching the AGNs and 
control galaxies to within a factor of 1.5 in PACS cover- 
age, in addition to M^^. This severely limits the number 
of matched comparison galaxies for some of the AGNs in 
our sample, but does not change our main conclusions. 
Any differences in the spatial distribution of the AGNs 
and the comparison galaxy sample on the sky does not 
lead to a lower fnd of the AGNs. 

Another possible source of bias arises from the use of 
MIR (MIPS 24/im) sources as priors during the photo- 
metric extraction of the PACS maps. In the small num- 
ber of cases where a PACS source was a blend of multiple 
MIPS sources, FIR fluxes were decomposed using the 24 
/im fluxes as weights. Since AGNs can be more MIR 
bright than inactive galaxies as a result of warm dust 
emission from the circumnuclear regions and torus, they 
could be assigned a stronger weight in blended sources, 
which can bias the FIR fluxes of AGN, as a population, 
to inaccurately high values. To test this, we repeated our 
analysis but only including AGNs that had no neighbor- 
ing 24 /im MIPS sources within a radius of 5". The 
results were unchanged, implying that extraction biases, 
if present, have a negligible effect on the FIR photome- 
try of X-ray selected AGNs in the CDFs. We conclude 
that the low fnd of AGNs is robust to systematic and 
data-related selection effects. 

Finally, we consider the effects of AGN-heated dust 
emission on the FIR detection rates of AGNs. We have 
shown that the majority of these AGNs have luminosities 
that are too low to boost a significant fraction of weakly 
star-forming hosts into the PACS detectable regime - 
if this were the explanation of the low quenching frac- 
tions, PACS-detected AGNs in the CDFs should gener- 
ally show very warm/blue FIR colors, which is not ob- 
served (Fig. 6). 



Are AGN in quenching galaxies? 



0.5 < z < 1.0 



1.0 < z < 1.5 



1.5 < z < 2.0 




-1 



2-1 1 

A log vL, [60 /j,m] (erg/s) 



2 -1 



Fig. 7. — A comparison of the Lso offset (ALso) of X-ray selected AGNs and inass-inatched inactive galaxies from the star-formation 
Mass Sequence. The statistical uncertainty in the distributions of the inactive galaxies are shown by the shading in the histograms - dark 
grey sections show the la uncertainty, due to the scatter in the population as well as small number statistics, while the light grey sections 
show 2(7. The dashed line at ALqo= corresponds to the center of the Mass Sequence. The AGNs show rather similar distributions to 
inactive galaxies. 



6. 



DISCUSSION: ARE AGN HOSTS A SPECIAL 

POPULATION OF GALAXIES? 



We have undertaken a detailed analysis of the FIR lu- 
minosities and detection rates of low and moderate lumi- 
nosity AGNs in the Chandra Deep Fields in the context 
of the population of inactive galaxies. Our key finding 
is the high FIR detection rates of AGNs, implying that 
they arc significantly less likely to be hosted by galaxies 
undergoing a special phase of star-formation quenching, 
contrary to the results of earlier studies. This differ- 
ence arises because, at the high masses of most AGN 
host galaxies, the UV-to-optical color is a poor proxy 
for SFR, which has complicated earlier studies of the SF 
properties of AGNs based on the use of color-mass or 
color-magnitude diagrams. 

We also find that low and moderate luminosity X-ray 
AGNs have similar distributions of FIR luminosities as 
inactive SF galaxies, consistent with the notion that host 
galaxies of such AGNs are drawn from a population of 
mostly normal galaxies. In a study of z ~ emission-line 
selected AGNs from the SDSS, Pasquah et al. (2005) find 
evidence for a small (0.2 dex) enhancement in mean FIR 
luminosities over a control sample of normal SF galaxies, 
matched to the AGNs by multiple structural and photo- 
metric criteria. Given the difference in the AGN selection 
method and matching scheme, we are unable to make a 
detailed comment on the difference between our results 
and those from Pasquali et al. (2005). However, a small 
offset of the magnitude uncovered in that work may still 



be consistent with our findings, given the substantial un- 
certainties on the ALgo distributions in Fig. 7 and the 
fact that the Figure excludes some weakly star forming 
galaxies. 

What could lead to the peculiar SF nature of AGN 
hosts? To answer this, we need to consider the rela- 
tive timescalcs of bright AGN activity (i.e, detectable in 
X-ray surveys of high rcdshift galaxies) and the modu- 
lation of SF on galaxy scales. It is highly unlikely that 
SMBHs accrete at a constantly high rate to maintain 
detectable AGN activity for a long period of time; typ- 
ical constraints on the lifetime of a bright QSO phase 
suggest 10^~^ years (Martini & Weinberg 2001; Hopkins 
ct al. 2005; Shankar et al. 2010). Even during a period 
of substantial accretion, the luminous output of an AGN 
may change quite significantly. High resolution simu- 
lations of SMBH accretion, both with and without the 
effects of feedback, suggest that the supply of gas to the 
accretion disk is not expected to flow in at a constant 
and steady pace (Hopkins & Quataert 2010; Novak et al. 
2011). Therefore, the X-ray bright AGN population is 
probably a rather transient population among galaxies 
in the field. 

On the other hand, SF on scales of the entire galaxy 
varies on timescales comparable to its dynamical time 
of around 10^ yr. In addition, the FIR luminosity, aris- 
ing mostly from dust heated by stars of a range of ages, 
is not a prompt measure of SFR, but can average over 
timescales of tens to hundreds of Myr, especially in sys- 
tems with fairly steady SF histories. 
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Fig. 8. — Histograms of the PACS non-detection fractions - the percentage of objects not detected in PEP+GH PACS maps - for 
1000 realizations of the mass-matched comparison sample of inactive galaxies in the corresponding redshift bins. The median value of the 
histograms is shown by the location of the vertical label 'Control'. The non-detection fraction of X-ray AGNs in the same redshift bin is 
shown as a thick arrow for comparison. AGNs have a significantly higher chance of being detected in the deep Herschel data, which implies, 
given the depth of the PEP-)-GH maps, that they preferentially avoid weaJfly star-forming, quenching or quiescent galaxies. 



In light of this, wc may consider two possible alterna- 
tives to explain our observations: 

A) Some evolutionary models suggest that a major 
fraction of moderate luminosity AGNs arise in the af- 
termath of a strong starburst, during which the SMBH 
also grows in a high accretion-rate QSO-likc phase (e.g., 
Ciotti & Ostrikcr 2007; Hopkins et al. 2008). A pop- 
ular form of these models suggests that gas-rich major 
mergers modulate both starbursts and QSOs. The end 
product is a quenching post-starburst galaxy hosting a 
low-to-modcratc luminosity AGN as the final amounts 
of gas fall intermittently into the SMBH. Within this 
scenario, the high PACS detection rates of AGN hosts 
may plausibly be ascribed to the gradual decline in dust 
heating by a post-starburst population, which serves as a 
fossil record of the co-evolutionary phase which occurred 
within the last few 100 Myrs. This picture would require 
a fine-tuning of the timescale between the peak of SF 
and the period over which low and moderate luminos- 
ity AGNs exist, in order to place the AGNs on the SF 
Mass Sequence with the same distribution of SSFRs as 
smoothly evolving inactive galaxies. It is also inconsis- 
tent with the observation that the [O II]-derived SFRs 
of moderate luminosity AGNs are comparable to normal 
star- forming galaxies (Silverman et al. 2009): emission 
line tracers are a more prompt measure of SF than the 
FIR and the starburst co-evolution scenario would there- 
fore predict a lower SF-powered line luminosity in AGNs 
than most comparable inactive galaxies. These observa- 
tions disfavor the post-starburst/major merger picture as 
the driver for most low and moderate luminosity AGN 
activity. 



B) A second view postulates that AGN hosts are sim- 
ply drawn from a smoothly evolving population of mas- 
sive star-forming galaxies. The connection here between 
star-formation and AGN activity is through the supply 
of cold gas needed for both, mediated through evolution- 
ary processes that modulate SF in galaxy disks and carry 
gas to the nucleus for eventual accretion onto the SMBH. 
The lack of any clear relationship between nuclear and 
SF luminosities in low and moderate luminosity AGNs 
(Shao et al. 2010; Mullaney et al. 2012b; Rosario et al. 
2012: Harrison et al. 2012a) argues against a prompt or 
synchronized connection between galaxy and black hole 
growth, favoring this scenario. The lower incidence of 
AGN among weakly star-forming or quiescent galaxies is 
a consequence of the depleted supply of cold gas in such 
systems; SMBH activity in such galax:ies is primarily in 
the form of a mechanically-dominated low accretion rate 
mode, such as that found in most radio-loud AGNs (Chu- 
razov et al. 2005; Trump et al. 2011; Best & Heckman 
2012). 

If cold gas is necessary for fueling X-ray bright AGN, 
then how do we reconcile this with the existence of such 

AGNs among genuinely quiescent galaxies, albeit at a 
low rate? This is indeed a puzzle, but is likely related 
to the fact that cold gas is not entirely absent in such 
systems. As many as 22% of nearby early-type galax- 
ies harbor substantial molecular gas (Young et al. 2011) 
and some show clear SF (Crocker et al. 2011). Dusty 
disks with very low SF efficiencies are commonly seen in 
the circumnuclear regions of massive elliptical galaxies, 
while filamentary dust is frequent on larger scales (Lauer 
et al. 1995; Ferrarese et al. 2006). Kauffmann & Heck- 
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man (2009) have suggested that SMBH accretion rates 
in quiescent galaxies are set by the supply of gas from 
stellar mass loss. Indeed, only a small quantity of this 
gas needs to intermittently reach the nucleus to fuel such 
low-luminosity AGNs. 

Existing surveys of large, complete samples of X-ray 
selected AGNs suggests the existence of a universal accre- 
tion rate distribution for SMBHs, which is independent 
of the stellar mass of the host (Aird et al. 2012; Bon- 
giorno et al. 2012). Our results can qualiiy these results 
by noting that the accretion rate distribution of AGN 
hosts do vary with SFR - galaxies with on-going star- 
formation are more likely to host a modestly accreting 
black hole than those that are quiescent. At lower stellar 
masses, both the gas supply and the SMBH mass drops, 
which is what possibly maintains the universal accretion 
rate distribution. However, the relative number of SF 
to quenched galaxies increases greatly with redshift — by 
2^2, 70% of massive galaxies are widely forming stars 
(Fontana et al. 2009). Concurrent to this is an increase in 
the gas fractions of galaxies (Daddi et al. 2010; Tacconi 
et al. 2010, 2012). Therefore, the accretion rate distri- 
bution must change with redshift, both in normalization 
and shape, if, indeed, the supply of gas is what limits the 
time-averaged accretion rate of SMBHs. The form of this 
change would be a shift to a higher characteristic break 
in the accretion rate distribution, as more SMBHs of a 
given mass accrete from the more abundant cold gas sup- 
ply at high redshifts. This is consistent with (MuUaney 
et al. 2012a) who show that average BH growth rates in 
mass-matched samples of galaxies indeed increase with 
redshift in line with their SFRs. 

As a final point, we consider the effects of galaxy dust 
properties, which could alter the interpretation of our re- 
sults. If AGN are preferentially in galaxies with larger 
dust masses, this would lead to redder UV-optical col- 
ors over galaxies with smaller dust masses or lower av- 
erage extinction. However, since dust is closely tied to 
gas, large dust masses would imply large gas reservoirs 
and, generally high SFRs, which only serves to under- 
line our findings. The difficulty of using U — V (or other 
similar colors) to understand SF properties in massive 
galaxies is compounded by the effects of dust extinction 
on these colors. For example, Cardamone et al. (2010) 
find that AGN hosts are common among galaxies where 
dust reddening has pushed colors from the star-forming 
Blue Cloud to the Green Valley. This qualitatively agrees 
with our findings that AGNs are predominantly in star- 
forming galaxies. Cardamone et al. (2010) also suggest 
that AGN hosts in the Red Sequence and in star-forming 
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galaxies have different accretion histories and feedback 
mechanisms. This may be true, but we maintain that, 
given the presence of molecular dusty gas even in gen- 
uinely quiescent galaxies, cold gas accretion may still be 
responsible for some of the AGN activity on the Red Se- 
quence. 

7. SUMMARY 

We characterize and study the star-forming properties 
of X-ray selected active galaxies in the Chandra Deep 
Fields North and South using deep Herschel/PACS pho- 
tometry at 70; 100 and 160 fim. Comparing the star- 
formation rate distributions of X-ray selected AGNs to 
those of inactive galaxies, we confirm that, after account- 
ing for stellar mass selection effects, AGN hosts lie on the 
Star- Formation Mass Sequence out to ~ 2, consistent 
with earlier studies (Mullaney et al. 2012b; Santini et al. 
2012b). However, we also find that AGNs are much more 
likely to be hosted by a star-forming galaxy than a qui- 
escent or quenching galaxy. This implies that AGNs are 
not preferentially found in galaxies that are undergoing 
a transformation from the Blue Cloud to the Red Se- 
quence. Instead, AGN hosts are drawn primarily from 
the population of normal massive star-forming galaxies. 
This may be interpreted in a scenario where radiatively 
efficient AGNs, such as those selected by X-ray surveys, 
require a supply of cold gas to sustain such phases of 
SMBH accretion. Combined with results which show 
the lack of any clear correlation between global SF and 
nuclear activity in low and moderate luminosity AGNs 
(Shao et al. 2010; MuUaney et al. 2012b; Santini et al. 
2012b; Rosario et al. 2012), our findings suggest that 
there may not exist any direct causal link between SMBH 
accretion and overall galaxy growth, but rather a more 
indirect relationship governed by the supply of gas, the 
fuel for both AGNs and star formation. 
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